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Abstract 

Climate change occurs when the Earth's energy budget changes due to natural or possibly 
anthropogenic forcings. These forcings cause the climate system to adjust resulting in a new climate 
state that is warmer or cooler than the original. The key question is how to detect and attribute climate 
change. The inference of infrared spectral signatures of climate change has been discussed in the 
literature for nearly 30 years. Pioneering work in the 1980s noted that distinct spectral signatures would 
be evident in changes in the infrared radiance emitted by the Earth and its atmosphere, and that these 
could be observed from orbiting satellites. Since then, a number of other studies have advanced the 
concepts of spectral signatures of climate change. Today the concept of using spectral signatures to 
identify and attribute atmospheric composition change is firmly accepted and is the foundation of the 
Climate Absolute Radiance and Refractivity Observatory (CLARREO) satellite mission being developed at 
NASA. In this work, we will present an overview of the current climate change detection concept using 
climate model calculations as surrogates for climate change. Any future research work improving the 
methodology to achieve this concept will be valuable to our society. 

Background 

The climate of the Earth's atmosphere is determined by the energy input from the Sun and the energy 
radiated from the Earth. Climate change occurs when the energy budget changes due to natural or 
possibly anthropogenic forcings. These forcings cause the climate system to adjust, sometimes 
engendering feedbacks that amplify or mitigate the forcings, resulting in a new climate state that is 
warmer or cooler than the original. The key question going forward is how to detect climate change. The 
measurement of infrared spectral signatures of climate change has been discussed in the literature 
relating to lower atmosphere (tropospheric) climate for nearly 30 years. Kiehl [1983] and Charlock 
[1984] noted that distinct spectral signatures would be evident in changes in the infrared radiance 
emitted by the Earth and its atmosphere, and that these could be observed from orbiting satellites. 

Since the pioneering work of Kiehl and Charlock, a number of other studies have advanced the concepts 
of spectral signatures of climate change [e.g., Goody et al., 1995; Slingo and Webb, 1997]. Harries et al. 
[2001] showed convincingly the increase in the greenhouse effect due to methane using measurements 
of infrared radiance spectra taken 27 years apart from satellite spectrometers. More recently, Huang 
and Ramaswamy [2009] have shown that distinct infrared spectral signatures associated with specific 
aspects of atmospheric structure and composition change (e.g., cloud height, tropospheric water vapor, 
C0 2 ) can be uniquely determined. These signatures can be used to determine changes in atmospheric 
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structure or composition can be monitored over decadal time scales. Today the concept of using 
spectral signatures to identify and attribute atmospheric composition change is the foundation of the 
Climate Absolute Radiance and Refractivity Observatory (CLARREO) satellite mission defined in the Earth 
Science Decadal Survey [National Research Council, 2007]. The CLARREO mission requires measurement 
of essentially the entire infrared emission spectrum of the Earth and its atmosphere, from 5 to 50 
micrometers (200 to 2000 wavenumbers). 

Current Challenge 

The CLARREO is currently in the pre-formulation phase of the mission. There are at present no 
measured, complete infrared spectra to simulate CLARREO measurements and to evaluate the 
fingerprinting techniques. To demonstrate the climate change detection and attribution concept, 
climate model simulation runs of the Canadian Centre for Climate Modeling and Analysis (CCCMA) 
model from the Cloud Feedback Model Intercomparison Project (CFMIP) are used as surrogates for 
measured spectra from which to assess climate change. Figure 1 shows schematically the procedure of 
processing climate model data for this study. The CCCMA was initiated with C0 2 of 280 ppmv (baseline 
C0 2 scenario for pre-industrial conditions), after it equilibrated, it was run for 20 years. Then the C0 2 
was gradually increased until it reached 580 ppmv (double C0 2 scenario), the model re-equilibrated and 
was run for additional 15 years. Outputs of the CCCMA model were 60-year time series of several 
responses. These responses were climate atmospheric effects (e.g., tropospheric water vapor, surface 
temperature, etc.) for all geographies that evolved over time. Only nine effects were determined to be 
the best-known major contributors to climate change, and therefore, they were considered in this study. 

The mean steady-state effects for the baseline and double C0 2 scenarios were obtained from years 1-20 
( Vib, V 2 B,..., V 9 B ) and years 46-60 ( Vid, V 2 D,..., V 9 D ), respectively. Typical post-processing of the 
mean steady-state atmospheric effects was to calculate a corresponding spectral radiance by a 
moderate-resolution radiative transfer model (MODTRAN). The spectral radiance calculation was done 
twice, one for baseline and the other for double C0 2 scenario. The spectral radiance difference between 
the baseline and double C0 2 scenarios was then obtained. It is the radiance difference that is used for 
detecting and attributing climate change between two climate states. This process of generating 
radiance difference is shown in Figure 2. 

Typically, the climate model resolves the globe into grid boxes and calculates atmospheric conditions 
for all grid boxes. The grid box size or resolution varies depending on climate models. Grid box of the 
CCCMA model has resolution of about 3.7 degree in latitude by 3.75 degree in longitude. Therefore, the 
model has 48 horizontal bands (180 degree/3.7 degree in latitude) and 96 vertical bands (360 
degree/3.75 degree in longitude), resulting in a total of 4608 grid boxes. Procedures of processing 
climate model data, generating their associated spectral radiances, and calculating radiance differences 
described above were repeated 4608 times, one for each grid box. 
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where 

V, (t) = a (4608 x t ) matrix for lower tropospheric cloud effects at time t for 4608 grid boxes 
V 2 (t) = a (4608 x t) matrix for middle tropospheric cloud effect at time t for 4608 grid boxes 
V 3 (t) = a (4608 x t) matrix for upper tropospheric cloud effect at time t for 4608 grid boxes 
V 4 (t) = a (4608 x t) matrix for C0 2 effect at time t for 4608 grid boxes 
V 5 (t) = a (4608 x t) matrix for stratospheric water vapor effect at time t for 4608 grid boxes 
V 6 (t) = a (4608 x t) matrix for tropospheric water vapor effect at time t for 4608 grid boxes 
V 7 (t) = a (4608 x t) matrix for stratospheric temperature effect at time t for 4608 grid boxes 
V 8 (t) = a (4608 x t) matrix for tropospheric temperature effect at time t for 4608 grid boxes 
V 9 (t) = a (4608 x t) matrix for surface temperature effect at time t for 4608 grid boxes 

V ib 9 V id = a (4608 x 1 ) vector for mean steady - state lower tropospheric cloud effect for baseline C0 2 and double C0 2 scenarios 
V 2b 9 V 2d = a (4608 x 1 ) vector for mean steady - state middle tropospheric cloud effect for baseline C0 2 and double C0 2 scenarios 
V 3b 5 V 3d = a (4608 x 1 ) vector for mean steady - state upper tropospheric cloud effect for baseline C0 2 and double C0 2 scenarios 
V 4b 9 V 4d = a (4608 x 1 ) vector for mean steady - state C0 2 effect for baseline C0 2 and double C0 2 scenarios 

V 5b 9 V 5d = a (4608 x 1 ) vector for mean steady - state stratopo spheric water vapor effect for baseline C0 2 and double C0 2 scenarios 

V 6b 9 V 6d = a (4608 x 1 ) vector for mean steady - state tropospheric water vapor effect for baseline C0 2 and double C0 2 scenarios 

V 7 b 9 V 7 d = a (4608 x 1 ) vector for mean steady - state stratopospheric temperature effect for baseline C0 2 and double C0 2 scenarios 
V 8b ? V 8d = a (4608 x 1) vector for mean steady - state tropospheric temperature effect for baseline C0 2 and double C0 2 scenarios 

V 9b 9 V 9d = a (4608 x 1) vector for mean steady - state surface temperature effect for baseline C0 2 and double C0 2 scenarios 

Figure 1 Procedure of processing the CCCMA climate model data 
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Figure 2 Procedure of generating radiance differences representing climate atmospheric effect change. 


As discussed above, we are looking for changes in climate caused by changes in nine different effects 
(although there may be many more to consider with real data). 

In order to quantify the contribution of each of the nine effects under study, nine one-effect-at-a-time 
suppression experiments were performed, assuming that the interaction effects are less likely to 
contribute to climate change. Table 1 lists all 9 suppression experiments. For example, Run 2 was an 
experiment that suppressed lower tropospheric cloud. The spectral radiance of the suppressed lower 
tropospheric cloud was calculated by inputting the mean steady-state of baseline C0 2 for lower 
tropospheric cloud while inputting the mean steady-states of double C0 2 for all other 8 effects into 
MODTRAN. This radiance was then subtracted from the double C0 2 radiance to get the signature of the 
contributing change due to the change in lower tropospheric cloud. 
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Table 1 Suppression experiments 


Run ID 

Description 

Inputs to MODTRAN execution 

Radiance difference obtained 

Available radiance 
difference array size 

Run 1 

a. Experiment with baseline 
C02 at 280 ppmv 

b. Experiment with double C02 
at 560 ppmv 

a. V1B5.-5V9B 

b. V1D5.-5V9D 

i?(VlD,...,V 9 D)-i?(VlB,...,V 9 B ) 

4608x1 

Run 2 

Suppression of lower 
tropospheric cloud effect V a 

V IB 5 V 2 D v ”9 V9D 

R 

VlD 

4608x1 

Run 3 

Suppression of middle 
tropospheric cloud effect V 2 

V ID , V 2B 5 V 3 D V M, V 9 D 

R 


4608x1 

Run 4 

Suppression of upper 
tropospheric cloud effect V 3 

V ID 5 V 2D 5 V 3 B 5 V 4 D V 9 D 

R 

V 3 D, V 3B ,...) 

4608x1 

Run 5 

Suppression of C02 effect V 4 

V ID v-5 V 3 D 5 V 4 B 5 V 5 D V 9 D 

R 

V4D, V4B,...) 

4608x1 

Run 6 

Suppression of stratospheric 
water vapor effect V 5 

V ID v-5 V 4 D 5 V 5 B 5 V 6D V 9 D 

R 


4608x1 

Run 7 

Suppression of tropospheric 
water vapor effect V 6 

V ID v, V 5 D 5 V 6B 5 V 7 D V 9 D 

R 

V 6 d, Vsb,...) 

4608x1 

Run 8 

Suppression of stratospheric 
temperature effect V 7 

V ID ,•••, V 6D 5 V 7 B 5 V 8D 5 V 9 D 

R 1 

(...,V7D,...)-4-V7B,...) 

4608x1 

Run 9 

Suppression of tropospheric 
temperature effect V 8 

V ID v-5 V 7 D 5 V 8B 5 V 9 D 

R 1 


4608x1 


5 




Run 10 

Suppression of surface 
temperature effect V 9 

V ID v .. 5 V 8 D 5 V9B 

R 


-R 



4608x1 

*(...) is s 

pectral radiance obtained from MODTRAN for a given climate mean steady-state effect inputs. 


is a spectral radiance 


vector for all mean steady-state effects, including effect i, from the Double C02 scenario. is a spectral radiance vector for the mean 

steady-state effect / from the Baseline C02 and all other mean steady-state effects from the Double C02 scenario. 
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CLARREO is interested in detecting and attributing climate change globally and zonally. Typical zone is 
defined as a 10-degree latitude band. Figure 3 gives an example of the spectral signatures for a 
particular zone. Figure 3a represents a measurement of the change in atmospheric radiance due to the 
doubling of carbon dioxide in the atmosphere. This is the radiance difference obtained from Run 1. 
Figures 3b to 3j show the spectral signatures associated with changes in lower tropospheric cloud, 
middle tropospheric cloud, upper tropospheric cloud, C0 2 , stratospheric water vapor, tropospheric 
water vapor, stratospheric temperature, tropospheric temperature, and surface temperature. These are 
the radiance differences obtained from Runs 2 to 10 respectively. They also aggregately explain the 
change in the infrared spectrum (Fig. 3a.) In practice, it is expected that the spectral changes (Fig. 3a) as 
a function of radiance vs. wavenumber (a reciprocal of wavelength) will be derived from the actual 
infrared satellite measurements. The spectral changes of all atmospheric effects (Fig. 3b-3j) will be 
derived either from the climate models under the control suppression experiments and from the 
existing satellites that were designed to measure specific atmospheric effects. Figure 4 shows the 
signature in each 10-degree latitude band varies across all 18 zones for two runs (Run 1 and Run 3). 




Wavenumber 


Wavenumber 


Wavenumber 


Wavenumber 


Wavenumber 


Figure 3 Example of the spectral signatures for a particular 10-degree latitude band (a) simulated 
change in infrared measurement, (b)-(j) spectral signatures due to changes in nine different atmospheric 
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parameters (e.g., water vapor, cloud height, temperature, etc.) aggregately explaining the change in the 
measure (a). 




(a) Simulated changes in infrared spectra due to doubling of carbon dioxide in Earth's atmosphere 
as a function in 10 degree latitude bands (18 total zones). 
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(b) Middle Tropospheric Cloud 


Figure 4 An example of (a) simulated changes in infrared spectral measurement (Run 1). (b) radiance 
changes associated with middle tropospheric cloud changes (Run 3). 


8 


There are two research objectives. The first objective is to identify globally and zonally which of the 
atmospheric effects change from the Baseline C0 2 to the Double C0 2 scenarios. The second is to 
quantify the magnitudes of all effects' changes. Achieving both objectives enables one to monitor 
atmospheric composition changes global and zonal over decadal time scales. The key challenge is how to 
achieve both objectives given many radiance difference spectral signatures are similar as shown in 
Figure 3. 

Future Research Work 

Future research focusing on exploring potential methodology is under way. The successful methodology 
will be able to detect globally and zonally which distinct infrared spectral signatures associated with 
specific atmospheric composition have changed. The accurate magnitude of change is also desired as 
this is an indication of which atmospheric component has changed. The current available data as 
described in Table 1 can be used to design and test potential methodologies. Future work to 
demonstrate this concept with the satellite spectral measurement is also planned. The proven successful 
methodology will be applied to the decadal time series of simulated-CLAREEO-like measurements to 
validate the utility of such measurements. 

Acknowledgement 

Harvard School of Engineering and Applied Science (Anderson research group) kindly provided the 
CCMA model output used in this study. 

References 

Charlock, T. P., (1984), C0 2 induced climatic change and spectral variations in the outgoing terrestrial 
infrared radiation, Tellus, 36B, 139-148. 1984 

Goody, R. M., R. Haskins, W. Abdou, and L. Chen, (1995), Detection of climate forcing using emission 
spectra. Remote Sensing and Earth Observation, 5, 22-32. 

Harries, J. E., et al. (2001), Increases in greenhouse forcing inferred from the outgoing longwave 
radiation spectra of the Earth in 1970 and 1997, Nature, 410(6826), 355-357, 
doi:10. 1038/35066553 

Huang, Yi, V. Ramaswamy, 2009: Evolution and Trend of the Outgoing Longwave Radiation Spectrum./ 
Climate, 22 , 4637-4651. doi: 10.1175/2009JCLI2874.1 

Kiehl, J. T., (1983), Satellite detection of effects due to increased atmospheric carbon dioxide, Science, 
222, 504-506. 

National Research Council, 2007, Earth Science and Applications from Space: National Imperatives for 
the Next Decade and Beyond, 456 pp., http://www.nap.edu/catalog.php7record id-11820 
Slingo, A., and M. Webb, (1997), The spectral signature of global warming, Q. J. R. Meteorol. Soc., 123, 
293-307. 


9 


